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I.  INTRODUCTION 


Waveguide  discontinuity  problems  are  very  common  to  microwave 
engineering.  Various  approaches  for  different  type  of  waveguide 
transitions  have  been  developed  including  the  conventional  lumped 
equivalent  circuits  and  field  theory  techniques ^ .  Following  the 
approach  of  Clarricoats*  which  is  a  field  theory  technique  led  to 
more  detailed  work  and  numerical  applications  by  Mielke^.  In  brief, 
this  technique  involves  matching  the  normal  mode  expansions  in  each 
waveguide  section  at  the  junction  between  sections.  A  large  matrix 
equation  is  obtained  which  is  then  solved  for  mode  amplitudes. 
Transmission  coefficients,  reflection  coefficients,  and  power  flow 
are  calculated  in  terms  of  these  amplitudes.  Mielke's  original  work 
only  allowed  an  abrupt  transitions  from  one  waveguide  to  another. 
Within  this  junction,  two  additional  restrictions  are  imposed.  One 
is  that  each  junction  can  have  only  one  input  and  one  output  waveguide 
Second  is  that  the  smaller  waveguide  must  lie  entirely  within  the 
boundaries  of  the  larger  waveguide.  Within  this  framework,  the 
following  coupling  coefficient  calculations  have  been  obtained 
analytically  and  implemented  in  a  computer  code: 

-  on-axis  symmetrical  modes  (circular  to  circular,  circular  to 
coaxial,  and  coaxial  to  coaxial  waveguide  transitions) 

-  on-axis  and  off-axis  for  all  modes  (rectangular  to  rectangular 
waveguide  transition) 

Subsequent  work  on  circular  to  rectangular  waveguide  transition  case 
has  been  implemented  numerically  by  Palevsky^. 


The  first  part  of  the  work  presented  is  to  introduce  a  tentative 
formulation  which  is  able  to  handle  the  transition  problem  of  the 
type  that  permits  a  conducting  wall  on  each  side  of  the  junction. 
However,  the  junction  is  still  restricted  to  be  an  abrupt  transition 
from  one  waveguide  to  another.  Also  each  junction  is  still  allowed 
to  have  only  one  input  and  one  output  waveguide.  The  modal  expansion 
technique  is  again  applied  here*.  This  formulation  is  not  yet 
implemented  and  its  derivation  is  described  in  Chapter  II. 

The  second  part  of  this  work  is  the  development  of  the  convergence 
criteria.  Accurate  results  depend  critically  upon  judicious  choice 
of  the  important  modes  via  symmetry  used  in  each  waveguide  section. 
Ideally,  an  infinite  number  of  modes  should  be  used.  In  practice, 
truncation  errors  caused  by  use  of  a  finite  number  of  modes  can  be 
reduced  if  these  modes  are  wisely  selected.  Criteria  for  such 
selections  are  derived. 

The  third  aspect  of  this  thesis  is  to  implement  the  evaluations 
of  the  coupling  coefficients  for  the  following  cases: 

-  on-axis  asymmetrical  modes  (circular  to  circular,  circular  to 
coaxial,  and  coaxial  to  coaxial  waveguide  transitions) 

-  off-axis  for  all  modes  (circular  to  circular,  circular  to 
coaxial,  and  coaxial  to  coaxial  transitions) 

The  numerical  implementation  of  the  above  cases  is  done  following  the 
technique  developed  for  the  circular  to  rectangular  waveguide  transition 


II.  FORMULATION 

The  transition  which  has  a  transverse  conducting  wall  on  each 
side  of  the  junction  is  investigated.  The  coaxial  to  coaxial  tran- 
sition  as  shown  in  Figure  2-lb  which  maintains  the  same  characteristic 
Impedance  is  a  common  transition  that  requires  this  extended  formulation. 
A  number  of  different  geometries  of  a  two-sided  boundary  transition 
is  illustrated  in  Figure  2-1. 

Adopt  the  definition  of  expressing  the  transverse  electromagnetic 
fields  in  the  expansion  of  an  infinite  number  of  vector  mode  functions 
found  in  Reference  3.  For  a  junction  which  has  an  input  guide  1  and 
an  output  guide  2,  the  fields  in  both  guides  are  matched  at  the 
junction. 

The  boundary  conditions  require  that: 
in  the  common  cross-sectional  area  (aperture) 

l't  (1)  -  \  (2)  (2-1) 

He  (1)  -  Ht  (2)  (2-2) 

and  on  the  left  side  of  the  conducting  wall 

Et  U)  -  0  (2-3) 

and  on  the  right  side  of  the  conducting  wall 

Et  (2)  -  0  (2-4) 


a  -  Circular  to  Circular  Transition 


Figure  2-1.  Two-Sided  Boundary  Transitions 
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Two  useful  pairs  of  equations  may  be  derived  from  equations  (2-1)  and 
(2-2)  by  manipulating  either  (2-3)  or  (2-4). 


First  scalar  multiply  equations  of  (2-1)  and  (2-3)  by 

N(l)  _  , 

£  e  (klf  1), 

k'  -  1 
1 


N(l) 


N(l) 


Et  (1)  *  V  e  (klt  1)  -  Et  (2)  •  T  e  (k1#  1)  (2-5) 

k  -  1  k  *  1 

1  in  the  aperture 


1 

-  NU)  -  ' 

Et  (l)  *  £  6  <kl*  °  "  0 

k'  -  1 
1 


on  the  left  side  of  the  conducting  wall 
Where  e  (k^,  1)  are  basis  vector  mode  functions  in  region  1  and 
N( 1 )  is  number  of  modes  used  in  the  expansion. 

N(2) 
k 


(2-6) 


__  I 


Secondly  post  vector  multiply  equation  (2-2)  by  £  e  (k2,  2), 

w'  «  1 


N<2>  _  ,  -  N(2)  _  , 

Ht  (1)  x  £  e  (k2.  2)  -  Ht  (2)  x  J  e  (k2,  2)  (2-7) 


k  -  1 
2 


k  ■  1 
2  in  the  aperture 


Where  e  (k2,  2)  are  basis  vector  mode  functions  in  region  2  and  N(2) 
is  number  of  modes  used  in  the  expansion. 

Integrating  (2-5)  over  the  aperture,  added  to  the  integration  of  (2-6) 
over  the  left  conducting  wall  gives: 


II- 


N  ( 1 ) 


> t  (1)  •  £  6  (kl»  l)  ds  "  lf  E t  (2)  ’  £  e  <k'l»  l)  ds 

k*  -1  1  k'  *  1 

cs ( 1 )  1  A( 1 , 2)  1  ( 2-8a) 

Integration  of  (2-7)  over  the  aperture  yields: 

ffL  N<2)  _  ,  1  _  rrL  n(2)  _  ,  1  _ 

JJ  |^Ht  (1)  x  Y,  e  (k2*  2)J  •  ds  -JJ  I  Ht  (2)  x  £  e  (k2,  2)J  •  ds 


N(l) 


A(l,2) 


A( 1 ,2) 


(2-9a) 


The  use  of  the  symmetric  manipulations  leads  to  the  following  dual 


equations  (2-8)  and  (2-9): 


Jit5' 


N(l) 

(i)  x  j 


A(l,2) 


kl  “  1 


e’  (kj,  1)]  .  ds"  “jfj^[Ht 


A(l,2) 


N(  1) 

(2)  x  7 

k '  ■  1 
K1  1 


e  (kj,  1)J  . 


ds 


ff N(2)  •  ff.  N(2)  _  . 

I  Et  (1)  .  £  e  (k2,  2)  ds  -  JJ  Et  (2)  .  £  e  (k2* 

A(l,2)  2  cs(2)  k2 


(2-8b) 

2)  ds 
(2-9b) 


Where  cs(l),  cs(2)  are  cross-sectional  area  of  guides  1  and  2  and 
A(l,2)  the  common  cross-sectional  area  at  the  junction  between  guides 
1  and  2.  Since  both  sets  of  equations  should  give  the  same  solutions 
in  the  limit  of  an  Infinite  number  of  modes,  any  linear  combination 
should  also  give  the  same  solutions.  How  these  two  sets  of  equations 
should  be  weighted  in  a  linear  combination  is  discussed  in  detail  in 
Appendix  A.  The  combined  equations  (2-8a  +  2-8b)  and  (2-9a  +  2-9b) 
form  the  matrix  equation  which  is  solved  for  the  forward  mode  amplitudes 
and  the  backward  mode  amplitudes.  The  size  of  the  matrix  will  increase 
as  the  number  of  modes  per  region  and  the  number  of  regions  Increase. 

To  set  up  the  large  matrix  equation,  the  coupling  terms  and  the 
propagation  terms  must  be  calculated.  The  one-sided  boundary  tran¬ 
sition  formulation  contains  coupling  terms  that  are  integrals  of  two 
modes,  one  from  each  side  of  the  junction,  over  the  aperture  (Overlap 
Integrals).  The  two-sided  boundary  transition  formulation  contains 
additional  coupling  terms  that  are  Integrals  of  two  modes,  both  from 
the  same  waveguide,  over  the  aperture  (Cross-Overlap  Integrals). 
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Whereas  Che  propagation  terms  are  the  same  for  both  transitions.  To 
include  this  type  of  transition  into  the  existing  system  of  equation 
of  the  one-sided  boundary  problem,  the  entries  of  the  matrix  structure 
as  well  as  the  entries  of  the  known  subvectors  have  to  be  modified. 

The  very  detailed  structure  of  the  matrix  equation  is  given  in  Appendix 
A. 


III.  NUMERICAL  CONVERGENCE  CRITERIA 


SUIManKVBWsi  mj*  w.  *■£*»" 


There  ere  always  solutions  to  the  matrix  equation.  However, 
unless  sufficient  modes  in  both  regions  preceding  and  following  any 
junction  are  considered,  the  solutions  will  be  inaccurate.  The 
accuracy  can  be  estimated,  in  principle,  by  examining  the  transverse 
fields  at  the  junction. 

This  is  done  by  adding  all  the  transverse  fields  in  region 
preceding  the  junction  to  obtain  the  total  field  and  substracting  the 
summation  of  all  the  transverse  fields  following  the  junction.  The 
result  will  indicate  how  close  the  fields  in  both  regions  are  matched 
at  the  junction.  However,  this  is  a  cumbersome  process.  Therefore, 
another  technique  is  developed  which  provides  a  way  of  estimating  the 
accuracy  of  the  solutions. 

Figure  3-1  illustrates  the  longitudinal  view  of  a  single  junction. 


Figure  3-1.  Longitudinal  View  of  a  Single  Junction. 


The  transverse  vector  aode  functions  of  guides  j  and  j+i  are  given 

by: 


e  (kj,  j)  for  guide  j 

e  (kj+i ,  j+l)  for  guide  j+l 


The  overlap  and  self-overlap  integrals  are  defined  by: 


0.1.  (kj ,  kj+1,  j)  •jj  S  (kj,  j)  .  e  (kJ+i,  j+l) 


(3-1) 


A(j,  j+l) 


S.O.I.  (kj,  j)  -  jje  (kj,  j)  .  e  (kj,  j)  ds 


(3-2) 


A(j,  j+l) 


S.O.I.  (kj+i  j)  -  jj e  (kj+| ,  j+l)  .  e  (kj+i,  J+l)  ds  (3-3) 


A(j,  j+l) 


First  consider  a  given  single  Incident  aode  in  larger  guide  j . 


Due  to  the  discontinuity,  an  infinite  nuaber  of  modes  in  guide  j+l 


are  excited.  This  incident  aode  over  the  saaller  cross  section  can 


be  represented  at  the  junction  by: 


N(j+i)  *  _  • 

e  (kj,  J)  -  lim  y  a  (kj+j,  j+l)  .  e  (kj+l,  j+l)  (3-4) 

N( j+l  )  ♦  •  k'  -1 

j+l 


Where  a  (kj+i,  j+l)  is  a  weighting  factor  of  each  excited  mode  In 
guide  J+l.  Substitution  of  (3-4)  into  (3-1)  yields: 


N( j+l ) 

0.1.  (k. ,  ki+l,  j)  -  lim  T  *  <kj+l»  J+D 

N(  j+l )  ♦  -  k'  L  -  1 

J+l 


.  jj  e  (kj+i ,  j+l)  .  e  (kj  +  i ,  j  +  l)  ds 
A(j,  j+l) 


(3-5) 
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Application  of  the  orthogonality  of  the  normal  mode  properties,  (3-5) 


becomes: 

0.1.  (kj,  kj+l,  j)  - 

/*  <kj+l 

(o 

Hence,  (3-4)  becomes: 

N(J+1) 

7  ( k 4 ,  j)  -  lim 

J  N( j+1 )  ♦  - 

Z 

kj+l  - 

Substitution  of  (3-7) 

1 

m 

o 

u 

e 

*  <kj+l.  J+D  if  kJ+l  “  kj+l 


If  k<+,  +  k4 


(3-6) 


(3-7) 


N(J+1)  ( 

S.O.l.  ( k  4 ,  J)  -  lie  V  jo.l.  (k1t  kj+j,  j) 

N(j+1)  +  -  kjTi  -  1  \ 


JJ  *  <kj+i»  *  •  (kJ+l»  J+1)  ds 


A(j,  j+1) 


(3-8) 


N(j+1)  (  )2 

S.O.l.  (k1t  j)  -  11m  T  (0.1.  (kj,  kj+1,  j)  (3-9) 

N(j+1)  ♦  -  kjT!  -  1  \  1 

Secondly  consider  a  single  mode  In  the  smaller  guide  j+1  Incident 

on  the  junction,  with  an  infinite  number  of  modes  are  excited  In 

guide  j.  This  incident  mode  over  the  large  cross  section  can  be 

represented  by: 


*  (kj+l»  J+D  *  li®  Z 

t 

N(j)  *  -  kj.i 


N(j) 

Z  b  (kj,  j)  .  e  (kj,  j) 


(3-10) 


Where  b  (kj,  j)  is  a  weighting  factor  of  each  reflected  mode  in  guide 
J.  Substitution  of  (3-10)  into  (3-1)  yields: 

N(j)  •  rr  _  - 

0.1.  (kj,  kJ+1,  j)  -  lim  £  b  (kJ»  j)  •//  e  (kJ’  J)  *  e  (kj»  J)  d* 

N(j)  ♦  -  kj  -  1  A(  J ,  j+1)  (3-11) 


(3-11) 


/b  (kj, 

>.I.  (kj,  kj+1,  j)  -  j 


j)  if  kj  -  kj 


if  k«  +  k* 


(3-12) 


Hence,  (3-10)  becomes: 


e  (k1+1,  j+1)  -  lim  Y  0.1.  (kj,  kj+1,  j)  .  e  (kj,  j)  (3-13) 

N(j)  -  -  k^l 


Substitution  of  (3-13)  Into  (3-3)  yields: 


e  (k.+1,  j  +  1)  •  e  (k. , , ,  j+1)  ds  -  lim 


cs(j+l) 


im  Y  (o.I.  (kj,  k  j+j ,  j)| 

N(j)  -  -  kjil  I  J  J  ' 

JJ  e  (kj,  j)  .  e  (kj,  j)  ds 


cs(j) 


(3-14) 


1  ■  lim 

N(j) 


N(j)  , 

Y  O.I.  (kj,  Kj+l,  J) 
♦  -  kj-1  1  ' 


(3-15) 


Equations  (3-9)  and  (3-15)  are  indicators  for  the  completeness 


of  the  total  representation  of  an  infinite  number  of  modes  accountable 


for  both  transverse  electromagnetic  fields  in  guide  j  and  j+1. 


For  the  first  case: 


N(j+1) 

£  {O.I.  (kj,  kj+1,  j)} 


kj+l  “  1 


S.o.l.  (kj,  j) 


(3-9) 


For  the  second  case: 


{O.I.  (kj,  kJ+I,  j)}  <1.0  (3-15) 


Therefore  the  convergence  can  be  verified  by  including  a  sufficient 
number  of  high  order  modes  in  both  regions  preceding  and  following 
the  junction  until  the  calculated  value  of  equations  (3-9)  and  (3-15) 
reaches  roughly  95Z  or  higher. 

The  knowledge  of  the  field  patterns  of  different  modes  for  various 
geometries  is  helpful  in  choosing  the  correct  modes  for  convergence. 

The  number  of  modes  required  for  convergence  varies  depending  on  the 
ratio  of  the  cross-sectional  areas  of  the  two  guides,  the  amount  of 
offset,  and  the  similarity  of  geometry.  For  example,  two  rectangular 
guides  of  only  slightly  different  dimensions  will  require  only  a  few 
modes  to  obtain  good  convergence.  However,  a  rectangular  to  circular 
junction  with  offset  axes  and  a  large  difference  in  cross-sectional 
area  will  take  many  modes  to  achieve  good  convergence. 
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IV.  NUMERICAL  EVALUATION  OF  OVERLAP  AND  SELF-OVERLAP  INTEGRALS 
FOR  CIRCULAR  AND  COAXIAL  GEOMETRIES 

For  a  one-sided  boundary  transition,  the  closed  form  solutions 
of  the  overlap  and  self-overlap  Integrals  for  the  symmetrical  modes 
on-axls  transition  are  available  in  Reference  4,  whereas  for  the 
asymmetrical  modes  on-axis  and  for  all  the  modes  off-axis  transition 
have  not  yet  been  obtained  analytically.  Therefore,  they  must  be 
evaluated  numerically.  The  overlap  and  self-overlap  integrals  for 
the  asymmetrical  modes  on-axis  transition  which  are  dependent  only  on 
the  radial  coordinate  are  evaluated  in  one  dimension,  while  the  over¬ 
lap  and  self-overlap  integrals  for  all  the  modes  off-axis  transition 
are  evaluated  in  two  dimensions. 

The  IMSL^  routines  DCADRE  and  DBLIN  are  used  for  one  and  two 
dimension  integrations  respectively.  DCADRE  is  a  routine  of  one 
dimension  numerical  Integration  of  a  function  using  cautious  adaptive 
Romberg  extrapolation;  where  DBLIN  is  a  two  dimensional  version  of 
DCADRE  applying  the  same  numerical  approach.  The  routines  require  as 
inputs  the  limits  of  integration  and  the  value  of  the  Integrand  as  a 
function  of  position.  The  remainder  of  this  chapter  describes  how 
these  two  items  are  calculated. 

Since  the  radial  dependent  part  of  the  mode  vector  functions  are 
expressed  in  terms  of  Bessel  functions,  an  explicit  evaluation  of  the 
functions  in  the  integrands  is  computationally  expensive.  Therefore, 
a  spline  routine  is  used  to  fit  the  radial  dependent  part  of  the  mode 
vector  functions,  and  spline  interpolation  is  used  to  evaluate  it  in 
the  integrand.  The  IMSL^  routines  ICSICU  and  ICSEVU  are  used  for  the 
fitting  and  interpolation.  Eight  steps  (radial  increment)  per  wave¬ 
length  plus  four  steps  at  each  end  point  are  used  in  the  spline  fit 
to  guarantee  a  smooth  interpolation. 
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v-v* 


b  -  Circular  -  Coaxial. 


guide  j 


guide  j  +  1 


c  -  Coaxial  -  Coaxial. 


Figure  4-1.  One-Sided  Boundary  On-Axis  Transition. 


The  integral  P  (kj,  k j+j ,  ^ )  is  evaluated  for  different  modes  coupling 
and  Che  results  are  listed  in  Table  4-1.  Obviously,  the  polarization 
factors  for  mode  couplings  of  the  different  azimuthal  mode  number  are 
equal  Co  zero.  The  remaining  task  is  to  numerically  Integrate  the 
radial  part  of  integrand.  The  limits  for  the  radial  Integral  are 
simply  the  extreme  boundaries  of  the  common  cross-sectional  area  at 
the  junction. 


Table  4-1.  Polarization  Factors  for  Various  Mode  Coupling  of  the 
Same  Azimuthal  Mode  Number  m  — 


Mode  Coupling 

Polarization 

TE  -  TE 

TE  -  TM 

— 

TM  -  TE 

TM  -  TM 

Parallel  - 
Parallel 

+  n  cos  (m.A9  ) 

-  it  sin  (m.A9  ) 

+  it  sin  (m.A6  ) 

+  x  cos  (in  .AS  ) 

Parallel  - 
Perpendicular 

-  *  sin  (m.A9  ) 

-  it  cos  (m.A9  ) 

s  cos  (m.A9) 

-  it  sin  (m.A0  ) 

Perpendicular  - 
Parallel 

+  x  sin  (m.A9  ) 

+  v  cos  (m.A9  ) 

-  w  cos  (m.A9  ) 

+  it  sin  (m.A8  ) 

Perpendicular  - 
Perpendicular 

+  it  cos  (m.AB  ) 

-  it  sin  (m.A9  ) 

+  w  sin  (m.AB  ) 

+  ft  cos  (m.A9  ) 

A0  is  an  angular  rotation  of  the  secondary  axes  relative  to  the  main 
frame  axes. 

B.  Overlap  Integral  of  an  Off-Axis  Transition 

The  off-axis  transition  involves  the  integration  of  a 
function  of  two  variables.  The  Integration  is  always  performed  in 
the  cartesian  coordinate  system  of  the  first  region  j.  The  relation¬ 
ship  between  the  two  coordinate  systems  is  illustrated  in  Figure  4-2. 
Ax,  Ay,  and  A0  are  the  offset  of  the  guide  j+1  relative  to  the  guide 
j,  and  A9  is  applied  after  Ax  and  Ay. 

16 


The  vector  node  functions  in  guides  j  and  j+1  can  be  written  in  terms 
of  their  cylindrical  coordinates  as: 


<kj,  j)  -  Er  (kj,  j,  r)  .  p  (kj,  J,  +) 


.  r«  + 


%  (kj.  J.  r>  •  j,  ♦)  . 

(■ 

«  (kj+j  •  j+1)  *  Er«  (kj+j,  j+1,  r')  .  p  (kj+i,  j+1,  ♦  ')  .  tq  + 

/v  I 

Ei|»'  (kj+i >  j+i •  r')  •  q  <kj+l,  j+1,  ♦’).  4»0 

^  A  A  f  A  • 

Where  rQ,  4>0  are  unit  vectors  along  r  and  directions;  and  r0,  <i>0  are 
unit  vectors  along  r*  and  V ’  directions.  The  expressions  p  (kj,  j,  *), 
q  (kj,  j,  *)  and  p  (kj+1,  j+l,  ♦’),  q  (kj+i,  j+1,  ♦  ')  are  polarization 
terms  of  guides  j  and  j+1  respectively. 

Using  the  following  coordinate  transformation. 


(4-4) 


x  -  r  .  cos  ♦  -  ♦  .  sin  + 
y  ■  r  .  sin  ♦  -  ♦  .  cos  t 


(4-5) 


the  transverse  fields  of  guide  j  evaluated  at  a  local  point  A  are 
given  in  the  cartesian  coordinate  system  (reference  frame)  as: 

Ex  (kj.  J)  -  (Er  (kj,  j,  r)  .  p  (kj,  j,  ♦)}  .  cos  - 

A  A 

{\  (kj,  j,  r>  .  q  (kj,  j,  *)}  I  .sin  \pA 


(4-6) 


Ey  (kj,  j)  -  { Er  (kj,  j,  r)  .  p  (kj,  j,  * )}  .  sin  U*A  + 

A  A 

{%  (kj,  j,  r)  .  q  (kj  j,  *)|  I  .  COS  l*A 


„v,,: 


Where  at  point  A,  rA  and  i|>A  are  given  by: 


12  2 
'  xa  ♦  yx 


(4-7) 


ip.  -  tan 


-1  *A 


At  the  sane  point  A,  the  coordinate  transformation  from  the  second  guide 
to  the  first  guide  is  given  by: 

xi  -  xA  -  Ax 

yi  -  yA  -  Ay 

r 2  r 

r'  -  -yxi  +  yi 


f  *2  *2 

xa  +  yA 


(4-6) 


4,1  - 


\p'  “  *Pl  “  48 


,-l  yA 


Applying  the  coordinate  transformation  of  (4-5),  the  transverse  fields 
of  guide  j+1  expressed  in  terms  of  the  reference  coordinate  system  are: 

Ex  (kj+1.  J+1)  “  lEr'  <kj+l»  J+l»  r>)  -P  (kj+l*  )+1»  •  cos  ^1  _ 

A  A 

{E^.  (kj+1  *  j+1,  r’>  *q  (kj+l»J+l»  <,,>l  •  sin^i 

A 

(4-9) 

Ey  (kj+1,  j+1)  I  -  {Er*  (kj+|,  j+1,  r’)  .  P  (kj+1,  j+1,  <'’)}  I  •  Sin  + 


{E^.  (kj+i,  J+1,  r’)  .  q  (kj+1,  j+1,  4'’)}  .  cos 


l. »  Nl^i  < 


The  lower  and  upper  limits  of  the  horizontal  end  points  are  found  to 

be: 


xlower 


+  Ax 


upper 


+  aj+1  +  Ax 


(4-11) 


While  the  limits  of  the  vertical  end  points  are  given  by: 


[2  2 

yiower<x)  -  4y  -  yj aj+1  -  (x  -  Ax) 


(4-12) 


12  2 

yupper<x>  ■  Ay  +  \aj+1  -  (x  -  Ax) 


B2.  Coaxial  to  Coaxial  Transition 


The  cross  section  view  of  the  junction  is  depicted  in 


Figure  4-4.  The  limits  of  the  horizontal  end  points  are: 


x lower 


+  Ax 


upper 


+  aj+1  +  Ax 


(4-13) 


The  common  cross-sectional  area  is  exactly  the  area  bounded  by  the 


inner  and  the  outer  walls  of  guide  j+1.  For  simplicity,  the  annular 


aperture  is  divided  into  two  regions  of  integration.  Figure  4-5 


shows  the  top  half  annular  ring,  where  the  dashed  line  indicates  the 


upper  bound  and  the  solid  line  indicates  the  lower  bound  of  the 


region  of  integration. 
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Figure  4-5.  Top  Half  Annular  Aperture. 


For  the  top  half: 


yiower(x) 


for  | x  -  A  x |  >  bj 


2  2 

Ay  +  \  bj+i  -  (x  -  Ax)  for  jx  -  Ax|  <  bj 


y upper 


12  2 
Ay  +  \  aj+l  “  (*  -  Ax) 


The  bottom  half  annular  aperture  is  illustrated  in  Figure  4-6. 


i  ’ j "j  •>  ’j  ’  ■''j  *j  •» 


'?  V  V  v  V 


L. 


Upper  Bound 


Ay  j _ 

0  1  Ax 


Region  of 
Integration 


Figure  4-6.  Bottoa  Half  Annular  Aperture. 


For  the  bottom  half: 


yiover(x)  •  Ay  - 


-  (x  -  Ax) 


(4-15) 


yupper^x^ 


for  |x  -  Ax|  J>  bj+i 


-  (x  -  Ax)  for  | x  -  Ax|  <  bj+j 


m 


•tt! 


'  * » v 

**:.*•>£ 


B3.  Circular  to  Coaxial  Transition 

Since  the  area  of  integration  is  bounded  by  the  inner  and 
outer  walls  of  coaxial  guide,  the  Units  for  both  horizontal  and 
vertical  end  points  are  the  same  as  the  previous  case.  The  same 
techniques  are  used  to  numerically  evaluate  the  self -overlap  Integrals 
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V.  EXPERIMENTAL  VERIFICATION 


The  coapuccr  code  developed  for  asymmetrical  modes  of  Che  on- 
axis  and  the  off-axis  transitions  has  been  cested  with  two  separate 
configurations.  The  two  dimension  numerical  integration  code  was 
tested  by  verifying  the  measured  return  loss  (reflection)  of  an 
off-axis  transition.  The  one  dimension  numerical  integration  code 
was  tested  by  comparing  the  experimental  results  of  the  return  loss 
of  a  centered  circular  aperture  of  finite  thickness  (obtained  from 
Reference  3)  to  the  numerical  results. 

The  ratio  of  the  cross-sectional  areas  of  the  two  guides  at  the 
junction  determines  the  preliminary  guess  of  the  ratio  of  the  number 
of  modes  considered  in  each  guide  at  that  junction.  Mode  selection 
for  limiting  the  choices  of  modes  in  the  expansion  to  those  which 
strongly  affect  the  wave  solution  was  done  by  observing  the  strength 
of  the  coupling  between  each  higher  order  mode  and  the  dominant  mode. 
The  convergence  criteria  was  satisfied  in  all  cases.  Violation  of 
this  criteria  yielded  inaccurate  results. 

A.  Three  Regions  Off-Axis  Circular  to  Circular  Transition 

An  illustration  of  the  three  regions  transition  is  depicted 
in  Figure  5-1.  Two  units  of  the  transition  were  fabricated.  The 
dimensions  for  each  unit  are  listed  in  Table  5-1.  Notice  that  the 
only  difference  between  the  two  units  is  that  the  length  of  the  middle 
region  which  is  of  11.0  mo  for  unit  A,  and  15.0  on  for  unit  B. 

A  block  diagram  of  the  experimental  arrangements  for  the 
reflection  measurement  is  illustrated  in  Figure  5-2.  The  equipment 
used  for  the  setup  is  listed  in  Table  5-2.  Calibration  for  100Z 
reflection  was  done  by  placing  the  shorting  plate  at  the  output  of 
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Table  5-1. 

aenslons  List  of  Unite  A  and 

B 

Unit 

lions  (ton) 

A 

B 

dl 

15.90 

15.90 

*1 

30.28 

30.28 

d2 

24.0 

24.0 

12 

11.0 

15.0 

d3 

15.90 

15.90 

^3 

30.28 

30.28 

4.064 


4.064 


Table  5-2.  Equipment  List  for  Reflection  Measurement 


HP  8690B  Sweep  Oscillator 

Wiltron  Model  560A  Scalar  Network  Analyzer 

HP  7035B  X-Y  Recorder 

HP  R422A  Crystal  Detector 

HP  Model  R532A  Waveoeter 

HP  Model  R752C  10  dB  Directional  Coupler 

TRG  Model  340-634  Mode  Filter 

TEjo  to  TEoi  Mode  Converter 

Shorting  Plate 


W  - 


vv 

'VV 


the  mode  filter.  Since  a  high  quality  matched  load  was  not  available, 
the  reflection  measurement  was  performed  by  leaving  the  output  end 
of  the  transition  open  to  air.  This  was  treated  as  an  18  dB  of 
return  loss  in  the  code  since  the  average  measured  return  loss  of  the 
unterminated  transition  was  approximately  18  dB  across  the  entire 
frequency  band. 

The  first  calculated  return  losses  which  are  obtained 
with  no  input  reflections  considered  are  not  in  good  agreement  with 
the  measured  results.  They  are  shown  in  Figure  5-3  for  Unit  A,  and 
in  Figure  5-4  for  Unit  B.  We  theorized  that  this  might  be  due  to  the 
mode  filter  not  absorbing  higher  order  reflected  modes  excited  at  the 
junction  between  the  mode  filter  and  mode  converter.  To  take  this 
into  account,  the  overall  length  of  the  first  region  preceding  the 
transition  was  adjusted  and  an  input  reflection  of  502  was  added  in 
the  computer  code  after  a  number  of  runs.  The  results  of  the  measured 
and  the  calculated  return  losses  of  the  two  units  are  shown  in  Figures 
5-5  and  5-6  respectively.  The  calculated  curves  for  both  units 
follow  the  measured  curves  closely  across  the  frequency  band.  At  the 
lowest  end  of  the  band  which  is  near  the  cutoff  frequency  of  the 
dominant  mode  ( TE ^ j  mode),  the  calculated  return  loss  is  higher  than 
the  measured  return  loss  for  Unit  A.  The  ripples  observed  across  the 
midband  for  Unit  A  and  uniformly  across  the  entire  band  for  Unit  B 
are  probably  caused  by  the  imperfection  of  the  mode  filter. 

B.  Centered  Circular  Aperture  of  Finite  Thickness  in  Circular 

Guide 

A  small  centered  circular  aperture  in  a  transverse  metallic 


plate  of  finite  thickness  (TEj^  mode  in  circular  guide)  in  circular 


ALCULATED 


MF.ASURF.D 


MEASURED 


guid«  may  be  Interpreted  as  a  shunt  inductance^.  The  cross  sectional 


view,  the  side  view,  and  the  equivalent  circuit  are  Illustrated  in 


Figure  5-7.  At  the  reference  plane  T  for  a  centered  circular  aperture, 


the  normalized  inductive  susceptance  is  given  by 


iWL  -  2.344 
4R  8.40  M 


(5-1) 


where 


(5-2) 


\2\  1/2 


(5-3) 


3.412Rj 


'  •  n 

the  reflection  coefficient  can  be  expressed  in  term  of  0  as: 


1  JL 

2  Y, 


©r 


(5-4) 


and  the  return  loss  in  dB  is, 


return  loss  »  -20  log^Q  (p ) 


(5-5) 


Measurements  of  the  return  loss  taken  at  X  *3.20  cm  in  a  circular 


guide  of  Ji.  in.  diameter  and  the  centered  aperture  metallic  plate  of 
16 


I _ in.  thickness  with  various  sizes  of  diameter  of  the  aperture  are 
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shown  in  Figure  5.5-3  in  Reference  3,  Section  5.5  page  243  to  246  in 


term  of  In  the  computer  code,  the  circular  aperture  of  finite 


thickness  is  modeled  as  a  circular  guide  with  smaller  diameter  than 


the  main  guide,  and  its  length  is  the  thickness  of  the  aperture. 


This  is  a  three  regions  on-axis  circular  to  circular  transition.  The 


fcV/'*.' 


K’V 


mm 


•  *!  \  >V 


diameter  of  23.812  mm  and  the  length  of  100.0  mm  are  the  dimensions 
of  the  first  and  the  last  regions.  The  length  of  the  middle  region 
is  0.794  mm,  and  the  diameter  varies  from  6.35  mm  to  23.019  mm.  The 
model  of  the  transition  is  depicted  in  Figure  5-8.  The  theoretical 
and  the  experimental  results  extrapolated  from  the  graphs  (Figure 
5.5-3)  in  Reference  3  are  compared  to  the  numerical  results  and 
listed  in  Table  5-3. 

Disregarding  the  results  at  the  point  size  aperture  and  at  the 
limit  of  the  size  of  the  aperture  approaching  the  size  of  the  main 
guide,  the  experimental  and  the  numerical  results  agree  very  well 
within  the  average  of  >2.5  dB  and  -0.04  dB.  Where  at  the  point  size 
aperture,  it  is  seen  that  the  numerical  result  is  about  0.009  dB  less 
than  the  experimental  result.  The  probable  cause  is  an  insufficient 
number  of  higher  order  modes  in  the  modal  expansion.  It  can  be 
explained  that  a  point  size  aperture  which  can  be  Interpreted  as  a 
Dirac  delta  function  has  a  spatial  Fourier  representation  of  an 
infinite  sum  of  harmonic  components  where  the  modal  expansion  technlqu 
starts  to  break  down.  Likewise,  at  the  limit  of  the  size  of  the 
aperture  approaching  the  size  of  the  guide,  the  numerical  result  is 
roughly  11  dB  higher  than  the  experimental  result.  It  suggested  that 
for  this  low  scale  of  VSWR's  reading,  the  measured  data  obtained  from 
Reference  3  was  inaccurate.  This  is  perhaps  due  to  the  insensitivity 
of  the  equipment  available  at  that  time.  Also  36  dB  and  47  dB  are 
equivalent  to  1.6Z  and  0.45X  reflection  coefficient  respectively. 

This  error  is  less  than  1.22  in  absolute  terms. 
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Table 

5-3.  Return  Loss 

of  a  Circular  Aoerature  in 

a  Circular 

Guide  of  15/16  in.  Diameter 

Return  Loss  (dB) 

d  (in) 

Theoretical 

Experimental 

Numerical 

Ratio  in  dB 

4 

0.0117 

0.0141 

0.0049 

-  0.0092 

5 

0.0506 

0.0180 

0.0260 

+  0.0080 

6 

0.1837 

0.1097 

0.1067 

-  0.0030 

7 

0.6189 

0.4299 

0.3644 

-  0.0655 

8 

2.0437 

1.1682 

1.0820 

-  0.0862 

9 

6.6540 

2.8036 

2.7454 

-  0.0582 

10 

5.5145 

5.8259 

+  0.3114 

11 

9.6205 

10.0786 

+  0.4581 

12 

14.6900 

15.5396 

+  0.8496 

13 

22.2399 

22.7566 

+  0.5167 

14 

28.6419 

33.3532 

+  4.7113 

14.5 

35.9187 

46.9358 

+11.0171 

^Numerical 

Note  that  Che  ratio  in  dB  is  defined  as  -20  logjo  _ 

p Experimental 
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VI.  CONCLUSION 


The  nodal  expansion  technique  was  found  to  be  capable  of  handling 
any  systen  of  uniform  waveguides  Interconnected  by  abrupt  junctions. 

The  computer  codes  developed  for  the  asymmetrical  modes  on-axis 
transition  and  all  modes  off-axis  transition  have  been  tested  and 
verified  through  two  separate  experimental  verifications.  The  results 
obtained  are  satisfactory.  Also,  the  numerical  convergence  criteria 
has  been  tested.  Violation  of  the  criteria  has  been  shown  of  giving 
inaccurate  solutions.  Formulation  of  the  two-sided  boundary  transition 
has  been  derived.  The  implementation  of  this  formulation  requires 
modification  of  the  structure  of  the  matrix  equation  of  the  existing 
code.  To  expand  the  capability  and  the  versatility  of  the  code, 
future  work  should  also  include  the  implementation  of  ridged  waveguide, 
dielectric  loaded  waveguide  and  multiport  junctions. 


ANN 


Appendix  A:  Derivation  of  the  Matrix  Equation 

A  very  detailed  derivation  of  aatrix  equations  of  a  two-sided 
boundary  junction  is  presented  in  this  appendix.  This  generalized 
transition  has  conducting  boundary  on  both  sides  of  the  junction. 

The  transition  discussed  in  Reference  4  was  restricted  to  only  one 
side  of  having  a  conducting  wall. 

A  single  junction  of  two-sided  boundary  is  illustrated  in  Figure 

A-l. 


Plane  of  junction  j 


|  guide  j+1 

guide  j 

\ 

' —  Right  Side  Conducting  Wall 

Left  Side  Conducting  Wall 


Figure  A-l.  Longitudinal  View  of  a  Two-Sided  Boundary  Off-Axis 
Transition. 


41 


y  /  y  y.  *• v-w.  '•  a  **.v.  .*  v. 


N  N  > 


m  t  wu  mi  in  rju  u  wiw  im.vuw.wui*  i  wwtvj  iw.iy.'KVHiw..vjvi.  yj  w.  *■.  *■_  v- 


Thc  guide  j  represents  the  input  guide  and  the  guide  j+1  is  the  output 
guide.  The  junction  is  viewed  as  one  transition  of  possibly  many 
transitions  preceding  the  guide  J  and  following  the  guide  j+1.  In 
general,  one  source  term  per  mode  is  allowed  anywhere  along  the  length 
of  each  section.  Figure  A-2  shows  the  forward  going  wave,  the  backward 
going  wave,  and  the  source  term  of  a  single  region  j. 


Loop  for  Exciting  Source 


B  +  B' 

s 


F 


Input  Plane 


+  F' 
s 


Output  Plane 


Figure  A-2.  A  Single  Region  j  Showing  the  Forward  Going  Wave, 
the  Backward  Going  Wave,  and  the  Source. 
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The  forward  node  amplitude,  F,  and  the  backward  mode  amplitude,  B', 
are  the  response  of  this  cavity  to  the  defined  source  terms  Fs  and 


Bs.  These  terms  F  and  B'  are  the  solutions  to  the  matrix  equation 
discussed  in  Chapter  II. 

The  fields  at  any  location  along  any  waveguide  region  can  be 
expressed  in  terms  of  the  sum  of  orthogonal  functions  weighed  by 
their  respective  mode  amplitudes  and  propagation  terms. 

For  a  single  junction  as  shown  in  Figure  A-l ,  the  total  transverse 
electromagnetic  fields  are  given, 


At  the  input  plane  of  guide  j+1; 


N(  j  +  1 ) 


Et  <j+0  -  2  F  <kj+l.  J+D  +  B  (kj+i.  J+O  +  Bs  (kj+l,  J  +  D  * 

»_  _  1 


Vi- 1 


e  (kj+i,  j+O 


(A-l) 


N( j+1 ) 


Ht  (j+1)  - 


I 

‘j.r 


Z(k 


j  +  1 


(F  <kj+l»  J  +  1>  -  B  <kj+l*  1  +  1>  - 


Bs  (kj+l*  J+1)  *  h  (kj+l-  J+1) 


(A-2) 


At  the  output  plane  of  guide  j: 


Et  (j) 


N(j) 

I 


V1 


F'  (k,,j)  +  B'  (kitj)  +  Fs  ( k j , j ) }  *  e  (k.,j)  (A-3) 


Ht  (J) 


N  ( j ) 

I 


kj“  1 


h  (k. , j )  (  ,  i 

z  ^71)  |F'  <kjJ)  -  B’  (kj.j)  +  Fs  (kj,j)| 


(A-4) 


Where  K  j  ,  Kj+j  are  mode  indicies  for  an  arbitrary  set  of  N(j),  N(j+1) 
modes. 
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Ac  Che  junction,  Che  transverse  fields  of  Che  Cwo  regions  must 
be  matched  ac  various  areas.  Let  cs(J)  and  cs(j+l)  be  Che  cross- 
secCional  area  of  regions  j  and  j+1  respectively.  And  A( J ,  j+1)  is 
Che  common  cross-sectional  area  at  the  junction  j  between  regions  j 
and  j+1.  In  the  common  cross-sectional  area: 

V 

£t  (j)  “  Et  (j+1)  (A-5) 

Ht  (j)  -  Ht  (j+1)  (A-6 ) 

On  the  left  side  of  the  conducting  wall: 

Et  (j)  -  0  ( A-7 ) 

On  the  right  side  of  the  conducting  wall: 

Et  (j+1)  -  0  (A-8) 

Two  possible  pairs  of  equations  may  be  derived  from  equations 
(A-5)  and  (A-6)  by  using  either  (A-7)  or  (A-8).  To  cast  the  equations 
in  a  more  useful  form,  the  following  manipulations  are  used.  First 
consider  the  use  of  (A-7).  Scalar  vector  multiplication  of  equations 


(A-5)  and  (A-7)  by 


N(j)  _  , 

X  e  (kj» 


j)  yields: 


Et  (j) 


N(  j  )  _  , 

X  e  (kj*j) 


V1 


N(j)  _  , 

Et  (j  +  1)  *  X  e  (kj.j) 
V  1 

(in  the  aperture) 


Et  (j) 


N(  j  )  _  , 

X  e  ( kj , j  ) 


V  1 


(A-9) 


( A- 1 0  ) 


(on  the  left  conducting  wall) 


N(  j+1 ) 


Pont  vector  multiplication  of  equation  (A-6)  by  J  e  (kj+j,  j+1) 


k  j+l“  1 


yields: 


N(  j+1 )  ,  N(j+1)  _  f 

Hr  (j)  x  T  e  (k1  +  1,  j+1)  -  Ht  (j  +  1)  x  J  7  (kj+l,  j  +  1)  (A-ll) 

.  »  .  .  •  _  . 


kj+l“  1 


kj+r  1 


(in  the  aperture) 


The  integration  of  equation  (A-9)  over  the  common  cross-sectional  area. 


added  to  the  integration  of  equation  (A-10)  over  the  left  conducting 


wall  gives: 


Cf  .  N(J)  _  •  rr+  N(j)  _  ' 

JJ  E,  (j)  *  I  e  (k1tj)  ds  -  JjEt  (j+1)  ‘  £  e  (k^j)  ds 

cs(J)  k>  1  A(  j ,  j+1)  k>  1 

j  J  (A-12) 


The  integration  of  equation  (A-ll)  over  the  common  cross-sectional 


area  yields: 


N(j  +  D  ,  N(j+1) 

JJ  [h*  (j)  x  £  e  (kj  +  1,  j+l)J  *  dl  -  JJ  [  Ht  (j+1)  x  £ 

A(j,  j+1)  kj+1-  1  A(j,  j+1)  k^,- 


k  ■  1 
j+i 


e  (kj+i»  j+1>]  *  os 
Substituting  ( A- 1 ) ,  (A-2 ) ,  (A— 3 ) ,  and  (A-4)  into  (A-12)  and  ( A— 1 3 ) 


( A— 13) 


gives : 

N(  j  )  N(  j  )  , 

I  I  {f*  (kjf  j)  +  B'  (kj.j)  +  Fs  (kj  ,  j  )J  * 


VI  V1 


II  e  (k1# 


cs(j) 


(kj , j)  •  e  (kj , j)  ds 


N(j)  N(j+1)  , 

I  I  (f  (kj+1,  j  +  1)  +  B  (kj+1,  j  +  1)  +  Bs  (kj  +  l,  j  +  1  )|  * 


kj-  1  kJ+1-  1 


II  e  (kj  * 


J)  *  e  (ki+. ,  j+1)  ds 


(A-12) 


A(j,  j+1) 


EM: 


m 


]&<,• 


l.'«; 


N(J)  N(j+1) 

Z  Z  0*I«  <kJ.  k1+l»J>  (F  (kJ.J> 


kj"  1  kJ*l"  1  Z  *kJ* 


-  •'  (kj»j)  +  F*  <kj,J)}  “ 


N(  j+1 )  N(j+1) 

T  Z  c*0,1‘  (kj+i*  kj+i»  j)  F  (kj+i»  j+1)  “  B  (kj+i»  j+1) 

.  •*  .  1  ■  1  ■-!  n  '  i.'n  I 


kJ  +  l"  1  kj+l”  1 


Bs  (kj+1.  j+1)}  (A-13a) 


Second  consider  the  use  of  (A-8).  Using  completely  the  symmetric 
manipulations,  the  following  dual  to  equations  (A~12a)  and  (A~13a) 
are  obtained. 


N(j)  N(j) 

Z  Z  C.O.I.  <kj.  kj  »  j>  |F  <kj.J> 


-  B  (kjt  j)  +  Fs  (kj,j) 


V  1  kj  -  1 


Z  (kj ,kj 


N(J)  N( j+1 ) 

Z  Z  0.1.  <kj»  kj  +  l»  J)  {F  <kj+l»  J+D  “  B  (kj+l*  J+1)  “ 
kj-  1  kj+1-  1  Z  (kj+1,  j  +  1) 


Bs  (kj+1,  J+1)}  (A-12b) 


N  ( J )  N(j  +  D  ,  , 

Z  Z  (F  <kj.j)  +  B  <kj.J>  +  Fs  (kj.j)|  0.1.  (kj,  kj+1,  j)  > 

kj“  1  k  j+1  1 

N(j+0  , 

Z  F  (kJ+l.  J  +  D  +  B  <kj+l»  j+1)  +  BS  <kj+l»  j+1 >J  (A-13b) 

kj+l“  1 

Both  pair  of  the  equations  (A-12a),  (A-13a),  and  (A-12b),  (A-13b)  can 


be  written  in  matrix  form  as: 


Ma  x  -  ba 


(A-14) 


Where  a  ■  a  or  b 


tS8& 


ff 

RW:; 


The  matrix  Ma  ,  vectors  x  and  ba  are  shovm  in  Figure  A-3.  Where 
the  subaatrices  and  subvectors  are  shown  in  Figures  A-4  through  A-16. 
For  the  case  of  single,  two,  and  J  junctions,  the  structures  of  the 
matrix  equations  are  shown  in  Figures  A-17  through  A-19. 

Both  these  sets  of  equations  should  have  the  same  solution  in 
the  limit  of  an  Infinite  number  of  modes.  The  similarity  of  these 
equations  to  the  single  sided  boundary  equations  derived  in  Reference 
4  is  demonstrated.  First  examine  the  equations  (A-12a)  and  (A-13a) 
in  the  limit  of  the  area  of  the  right  wall  going  to  zero. 

Since, 


limit 

area  of  right  wall  * 


.// 


e  Oc<+i  ,  j+D  *  e  (k.+1 ,  j+1)  ds  - 


A(j,  j+1) 

Jfe  (kj+1,  j+1)  *  e  (kJ+1,  j+1)  ds 

cs(j+l) 


j  1  if  kj+1  -  kj+1 
1  0  ^  kj+1  +  kj+i 


(A-15) 


Then,  equations  (A-12a)  and  (A-13a)  become  a  pair  of  equations 
derived  for  larger-to-smaller  transition  defined  in  Reference  4. 
Similarly,  it  can  be  shown  that  in  the  limit  of  the  area  of  the  left 
wall  going  to  zero,  the  equations  (A-12b)  and  (A-13b)  are  equivalent 
to  another  pair  of  equations  derived  for  sraaller-to-larger  transition 
in  Reference  4. 

Since  both  sets  of  equations  should  yield  the  same  result,  any 
linear  combination  should  also  yield  the  same  result.  It  is  expected 
that  the  solutions  to  the  equations  (A-12a)  and  (A-13a)  should  converge 


Equation  (A12  a)  in  Matrix  Form: 


(a  -  a  or  b) 


o’o  7  .  (i+r 


F  (k  ,J) 


B*  (k  ,j) 


2N(j) 


F  (kj+1.  j+D 


B'  <^j+l.  j+1) 


2N ( j+1 ) 


Figure  A-16.  Solution  Subvectors. 
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faster  to  the  correct  solutions  as  the  number  of  modes  per  region  is 
Increased  for  Junctions  that  have  most  of  the  conducting  wall  on  the 
left  side.  Similarly  one  would  expect  the  solutions  to  equations 
(A-12b)  and  (A-13b)  to  converge  more  rapidly  to  the  right  answers  for 
junctions  with  most  of  the  conducting  wall  on  the  right  side  of  the 
boundary.  Therefore,  it  would  seem  reasonable  to  make  a  linear 
combination  of  equations  (A-12a),  (A-13a)  and  (A-12b),  (A-13b)  where 
each  term  is  weighted  by  the  relative  area  of  the  conducting  wall  on 
each  side. 

Thus, 


u>l  *  Ma  x  ■  o>£  *  bj 


u>. 


*  Mjj  x  •  ur 


(A-16) 


bb 


Where  w*  and  u)r  are  defined  by: 


area  of  left  wall 


area  of  left  wall  +  area  of  right  wall 


area  of  right  wall 


area  of  left  wall  +  area  of  right  wall 


( A— 17) 


After  combining,  it  yields: 

jwj  *  Ma  ♦  <*>r  *  M{j}  x  ■  ui£  *  ba  +  u)r  *  bj,  (A-18) 
Where  the  submatrices  and  subvectors  of  (A-18)  are  shown  in  Figures 
A-20  through  A-26  for  a  single  junction. 

The  backward  mode  reflection  coefficient  for  the  first  region 
that  generates  a  forward  going  wave  In  the  first  region  (not  a  source 
term)  is  normally  defined  for  a  given  problem.  Similarly,  the  forward 
mode  reflection  coefficient  for  the  last  region  that  generates  a 
backward  going  wave  in  that  last  section  is  also  a  known  quantity. 


z  <k| 


(Z*lx)  z 


This  Is  equivalent  tc  defining  the  input  and  output  termination  of 
the  problem.  Therefore,  the  structure  of  the  matrix  equation  for  J 
junctions  (J+l  regions)  is  modified  as  shown  in  Figure  A-27 ;  where 
the  modified  submatrices  and  solution  subvectors  are  shown  in  Figures 
A-28  through  A-32. 


< k i , 1)  •  Cf  (kj.l) 


( k i ,  1 )  .  VRI  (kj.l)  +  Uf 


Appendix  B:  Derive  ion  of  the  Self-Overlap  Integral*  for  Symmetrical 

Modes  of  Circular  end  Coaxial  Geometries  On-Axis  Transit  ion 

The  criteria  for  the  convergence  check  established  in  Chapter  III 
requires  a  knowledge  of  the  values  of  the  self-overlap  Integrals  of 
every  mode  accounted  for  at  the  Junction.  The  area  of  Integration  is 
always  the  common  cross-sectional  aperture.  The  self-overlap  Integral 


for  any  node  in  the  sculler  guide  is  trivially  unity.  Therefore,  only 
the  self-overlap  Integral  for  a  node  in  the  larger  guide  over  the 
common  aperture  is  derived. 

The  S.O.I.  for  any  mode  kj  in  the  larger  cross  section  guide  j 
is  defined  by 

S.O.I.  (kj.J)  -  jfj^e  (kj,  j)  *  e  (kjt  J)  ds  ( B—  1 ) 

A(J,j+l) 

The  closed  form  solutions  of  the  self-overlap  Integrals  for  the 
symsutrlcal  (m-0)  modes  on-axis  transition  are  derived  below  for 
various  waveguide  geometries. 

B.  1 .  Circular  to  Circular  Transition 

The  longitudinal  view  of  the  transition  is  depicted  in 
Figure  B-l.  From  the  Appendix  C,  the  TE  mode  vector  function  in 
guide  j  with  kj  representing  m  •  0  and  arbitrary  n  is  given  by: 


-TE  , 

e  ( k  j  ,  j  )  - 


n-  Jo  <*on>  *  V 
*  Sj  .  J0  (*on) 


(B-2 ) 


$ 


Where  n  ■  1,  2,  3,  ...  and  xon  is  the  nth  zero  of  JQ 
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$ 


Substitution  of  (B-2)  intlo  (B-l)  yields: 

r2'  f‘]*‘  /  n~  Jo  •  =j)l 

S.O.I.  (k4,J)  -  /  J  V”“  • - - T^T  •  rdr  d*  (B-3) 

J  JQ  ■'O  *  *  *J  •  Jo  <*on> 

Application  of  the  derivative  formula  of  Bessel  function: 


J„  (x)  -  -  Ji  (x) 

And  with  the  use  of  the  integrals  of  Bessel  functions  formula  (4)  on 
page  255  in  Reference  2,  (B-3)  becomes 


1  v  on  ‘j  r  0  l  on  »j  ) 

•  J2  V  °n  aj  , 

j2  (*  ) 
o  ^  on/ 

(B-4  ) 

For  TM  mode  with  kj  representing  m*0  and  arbitrary  n,  the  vector  mode 
function  is: 


-TM 

e  (kj 


/  1  Jo  (*on  * 

.D-W-  •  — v  .  H  ■% 


aj  .  Jl  (xon) 

Where  xon  is  the  nth  zero  of  J0 
Substitution  of  (B-5)  into  (B-i)  gives: 

('  f**'  \  Jo  ("on  •  5j) 

S.O.U  .  j  J  )  ,,  uj 

0  0 

Evaluating  the  Integration  yields: 


(B-5) 


r  dr  d<’  (B-6) 


/Vl\2  J1  (*on  ap)  -  Jo  (xon  ip).  J2  (*on  7p) 
1.0.1.  (kj,J)  "I - 1  .  -  - 
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B.2.  Coaxial  to  Coaxial  Transition 


Tha  type  of  transition  (one-sided  boundary)  as  shown  in 


Figure  B-2  is  valid  for  the  derivations  below.  For  the  case  of  kj 


representing  n*0  and  arbitrary  n,  the  TC  node  vector  function  is: 


_TE 

e  (kj,J)  - 


«  ■-  —  -  r  /  i  r  \  i  i  •  /  •  *r  \  •  • 

lon  V  *  Jo  ^xon  bjj.  No  ^xon^  “  No  yxon  bjj  .  ^o  ^xoi 


.2  ,• 

Jo  ^xon) 


j’2.  (/  aA 

o  xon  — 

\  bi) 


(B-8) 


Substituting  (B-8)  into  (B-l)  and  collecting  terns  yields: 


S.O.I.  (kj,J)  ■ 


y2  L-  *J 

Jo  l  xon  — 


No  (xon) 


Evaluating  the  Integral  yields: 


\  /  /  v 

-  N0  ^xon  —  ^  •  J0  ^xon^|  .  r  dr 


(B-9) 


S.0.1.  (k, ,J )  - 


1 2  , 

Jo  ^xon- 


f-2  U' 

o  I  xon  •  bj  1 


on  •  «i 


■)Vf, 

. 


'  /  ’  *J+1\ 

*  N1  <xon>  "  N1  (^xon  ‘  — j  *  J 


1  v  *on • 


^  X/in  )  /  -<N 


/  '  aJ+l\ 

[Xo*  ‘  —)•  Jl  (xor 


P 

jC 


*  ■*  1 
VVV' 

■'w 

'‘C<VV‘ 

•  «_  \v 


.s  •  /• 

v!  .■_*«*  >■/ v.v.  .•>}»•  .v.y.v.  ^ ■■.v.v,.' ^ .  /  .-v.  '"Sly-yl?  ■’ * 


aj+l 

bJ 


^*oo  •  •  * 1  <xoo>)  •  |*2  ^xon  bj~)  ‘  J*  (X°n>  ’  ^  (X°n  * 


•  *1  (*™>)  ~(^y  *  bJ+)  *  {Jl  (*°n  ’~rf)  * Nl  (x°n)  ' 

2 

/  *  bj^i\  i  /  *  bj+i  \  ’ 

*  1  y*on  •  — — J  •  -^1  ^xon^  “  |*o  ^xon  •  — j  •  ^1  (xon) 


( *on  • 


/  '  bj-i\ 
V  ■  “/  ' 


bj*l\  '  /  '  bj+l  \ 

r —  I  •  *1  (xon)  •  *2  [xon  •  — r - ]  *  ^  1  ^xon)  " 


*1  (xon) 


(B-10) 


For  TM  node  wlch  k.j  representing  a-O  and  arbitrary  n,  the  vector  mode 


function  is: 


-  TM  xon  VT  Jo 

*  ‘“j-11  *tf  l  - 


o  (xon  bj )  No  ^xon^  “  No  (xon  bJ ) 


(xon) 


Jo  ^*on  .  -5J-^ 


Substitution  (B-ll)  into  (B-l)  gives: 


2n  aj+l 

i.o.i.  (kj.j) . j  j  .  jy  . 

Jo  Jbj  +  1  bj  >/  2 

Jo  (xon  ^j)  •  "o  <xoa>  *  No  (xon  ’  *j)  •  Jo  <>on'  I 


J n 


J°  i'm  '  ^") 


•  Jo  <xo„> 


(B-ll) 


' ••••>..  v. 


Evaluating  the  integration  yields: 


S.O.I.  (kjtj) 


JQ  ^xon^ 


lon  • 


“?)  [(Nl  (*“  -  f1)  • 


Jo  (xon)  ” 


(J°  X)  ) 

Jo  (xon)  "  Jl  ^xon  •  ^  •  No  (xon^  j  “  j^o  ^xon  •  — — — ^  •  Jo  ^xon^  " 

Jo  ^xon  * - ^  *  No  (xon)  |  *  ^2  ^xon  •  — — ^  •  Jo  ^xon^  “  J2  ^xon  •  — — — ^ 

No  (xon>  jj  -  ^xon  |jNl  (*on  •-J— )*  J°  (Xon>  “  Jl  (X°n 

Nq  (*on)  |  “  |no  ^xon  •  — - — Jo  ^xon)  “  Jo  ^xon  •  — - — ^  •  N0  (xon)  |  • 

|n2  ^xon  *  — *  Jq  (*on)  "  J2  ^  xon  •  — - — ^  •  N0  (xon^  j^j  (B-13) 


Nq  (*on) 


For  TEM  mode  with  kj  representing  m»0  and  n-0,  the  vector  mode  function 


_  TEM 

e  (kj.j)  - 


(B-14) 


r  .  \l2n  In 


Equation  (B-l)  becomes: 


S.0.1.  (kjfj) 


■17 


v  .  2*  In 


r  dr  cty 


(B-l 5) 


ym; 


V 


',5 : 


i’V'v 


*0 


?A 

v, 

» 


Mi 


5» 


Evaluating  the  integration  yields: 


in 


S.O.I.  (kj.j) 


m 


in 


aj 


(B-16) 


B.3.  Circular  to  Coaxial  Transition 

The  derivation  below  is  valid  for  a  smaller  coaxial  guide 
totally  embedded  inside  a  larger  circular  waveguide.  Such  a  transition 
is  depicted  in  Figure  B-3. 

With  kj  representing  m«0  and  arbitrary  n,  the  vector  mode 
functions  for  a  circular  guide  are  given  in  (B-2)  and  (B-5)  above. 
Performing  the  Integration  over  the  annular  aperture  yields  expressions 
of  S.0.1.  For  TE  mode: 

2 


S.O.I.  (kj. 


4  <-;»>  lv  1 1  K  'i 1  v  aJ ' 


'  aJ+l  \1 

/ v>  \2  r, 

2  /  *  bj+l\ 

T  (  ' 

v  °n  *  J  r 

1  .  J  *  J 

1  l  xon - J 

-  Jo  l  xon - I* 

\  /J 

\  hj  /  L 

\  *j  / 

V  a  j  ' 

(-  •?)]) 


( B— 1 7 ) 


Where  xon  is  the  nth  zero  of  JQ 
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For  TM  node: 


S.O.I. 


'  ri: '  Kt)  '  I'!  ('“ '  v) '  (-  •  ?, 

(.00  •  JJi)]  - 

(Xon  ‘-I1)]} 


/ bJ+l\ 

ft2  /  bj^l\ 

f  /  VA 

l— -  • 

I  J1  (  xon  *  -  -  ) 

**  ^o  (  xon  •  - 1 

\  bj  / 

L  V  *j  / 

\  / 

(B-18) 


Where  xon  is  Che  nth  zero  of  J0. 


All  Che  above  incegrals  are  evaluaced  vlch  Che  aid  of  Che 
jincegrals  of  Bessel  funecions  formula  (4)  on  page  255  in  Reference  2. 


I 


» 
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Appendix  C:  Vector  Mode  Functions  for  TE.  TM.  and  TEM  Mode  for 
Rectangular.  Circular,  end  Coexlal  Geometries 

The  tranaverae  fields  for  the  rectangular,  circular,  and  coaxial 
waveguides  are  derived  in  Reference  3.  The  formulas  for  the  TE,  TM, 
and  TEM  modes  are  given  below.  The  origin  of  the  coordinates  have 
been  moved  to  the  center  of  the  guide  from  the  lower  left  hand  corner 
for  the  rectangular  waveguide  modes.  The  vector  mode  functions  are 
represented  by 

_IM 

ejm  =  TM  mode  vector  function 

_TE 

emn  -  TE  mode jvector  function 

_IEM 

eQQ  =  TEM  mode  vector  function  (coaxial 
guide  only) 

For  a  rectangular  waveguide,  the  cross  section  view  is  depicted 
in  Figure  C-l. 


y 


Figure  C-l.  Cross  Section  View  of  a  Rectangular  Waveguide. 


•*  »■* 


The  TM,  and  TE  vector  node  functions  for  a  rectangular  waveguide 


are  given  by 


2.  f|Vr  air  a 

»  IL9~  co®  a”  +  2^  •  8in 


nir  b  | 

F"  <y  +  7>J  •  *0 


0  b  9  a 
■  T  +  n*  x 


[nir  air  a  nir  bl  ^  I 

b  sin  a  (x  ♦  f)  .  cos  b  <y  +  2J  •  y0f  (C-l) 


ea  % 


emn  “ 


*  ( rhtf  air  a  nir  bl 

2  o - Ta  llF  cos  *  (*  +  2>  •  *in  b  <y  +  V\  •  *0 

»  a  +  n  F 

[bit  bit  a  nfr  bl  A  | 

■j-  sin  —  (x  +  -j)  •  cos  -g"  (y  +  j)J  .  y0[  (C-2) 


Where  xQ,  yQ  are  unit  vectors  along  x  and  y  -  axis  respectively 
a,  n  ■  0,  1,  2,  3,  ,  node  a  ■  n  ■  0  excluded 


1  if  a  -  0 


2  if  a  +  0 


Figure  C-2  illustrates  the  cross  section  view  of  a  circular  waveguide 


Figure  C-2.  Cross  Section  View  of  a  Circular  Waveguide. 


_TM  - 

« 


The  TM  node  vector  function  for  a  circular  waveguide  is: 

'iSMH  <'-'3 '[“]"]•  *•* 

[rfc ' J-  <-  '  9 '  [ -]  "] ' 


(03) 


The  TE  mode  vector  function  for  a  circular  waveguide  is: 


Bi  m  0|  1 1  2 1  ••• 
n  m  1|  2|  3)  ••• 

xmn  “  nth  non-vanishing  root  of  the  mth  -  order  Bessel 


function  Jm  (xnm)  -  0 

« 

xmn  “  nth  non-vanishing  root  of  the  derivative  of  the  mth  - 

•  I 

order  Bessel  function  JB  (x^)  -  0 

(1  if  m  -  0 
2  if  i  i*  0 


The  TM  mode  vector  function  for  a  coaxial  waveguide  is: 

JTH  V*£«  If X*1  C**1  *  T>)  ^x*1^  *  **m  (*mn  •  T>)  Jm  (*mn) 


(x__ ) 


M m_jw_  _  2\ 
Vm  <C  x^)  ) 


’]-*] 


The  TE  mode  vector  function  for  a  coaxial  waveguide  is 
m  Jm  {*mn  •  Nm  ^*mn^  “  Nm  fxmn  • 


\u«(cx«,))  i  \c  **/;  i  \xmn/  J/ 

'  *  /  *  r\  *  '  '  /  '  Jf  \  *  ' 

:mn  \xmn  •  b  /  Nm  ^xmn^  “  \xmn  •  b  / -*m  ^xmn^ 


■  -sin  -|  “I  ^  xmn  (xmn  •  b)  Nm  ^xmn^  “  (xmr 

ro  +  ~B - 

LcosJ  (j4^r  f-f-vt 

\  (  ^m  xmnV  \C  xmn/ 


P  -1  •  ;-l 


(C-6) 


^fe:45- 

M 


& 
v 


m 


V  ?/*■  v-  V-  v1.  V- 


And  the  TEM  mode  vector  function  is: 

—TEM  1 

•oo  - - - - —  •  'o  <M> 

r  .  \  2»  in  (£) 

A  A 

Where  r0,  \pQ  are  unit  vectors  along  r  and  directions 
m  m  0|  If  2f  3| 

n  -  1,  2,  3,  ... 


nth  non-vanishing  root  of  the  mth  -  order  Bessel- 


Neumann  combination. 


nth  non-vanishing  root  of  the  derivative  of  the  mth  - 
order  Bessel-Neumann  combination 


1  if  m  ■  0 


2  if  m  j*  0 
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A(J,  j+1) 
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B  (k,j) 

B'  (k, j ) 

(k, j ) 
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s 
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A( j ,  j+1) 


I 

e  (k j , j  >  ds 


f 

C.O.I.  (kjfj  t  kj+j ,  j) 


Cross-overlaop  integral  of  mode  kj+j  to  mode 
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common  cross-sectional  area  of  regions  j  and 
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i 

e  (kj+1,  j+1)  ds 
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electric  field  magnitude  of  mode  kj  in 
region  j 

The  radial  dependent  part  of  the  azimuthal 
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